Effects of mechanical strain on thermal denaturation of DNA 
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A class of simple statistical mechanical models for DNA melting, first proposed by Poland 
and Scheraga, have been demonstrated to exhibit first or second order thermodynamic singularity, 
notwithstanding the intrinsic one-dimensional nature of the problem. In the present paper we show 
that the inclusion of twist elastic energy in the Poland-Scheraga models either leads to suppression 
of the thermodynamic singularity or to a weak, third order singularity. 

d ■ 87.14.Gg, 05.70.Fh, 63.70.+h, 64.10.+h 

The thermodynamic properties of DNA near the melting, or denaturation, point has become the focus of intense 
interest in the theoretical physics community A class of simple models, proposed in their original form by 

Poland and Scheraga (PS) [Q exhibits a thermodynamic singularity, notwithstanding the essentially one-dimensional 
nature of the DNA molecule. In these PS models the single-stranded DNA sections inside denaturation bubbles are 
treated as excluded-volume polymers performing a random walk. This produces an effective long-range interaction 
that leads to a thermodynamic singularity in this one-dimensional system. Depending on whether excluded volume 
effects between single and duplexed strands are or are not taken into account, the thermodynamic singularity can 
i_S , have the character of either a first Q or a second-order || phase transition. We denote those two alternatives as Case 
I and Case II, respectively. 

Experimental studies of DNA melting have not yet provided any clear support for the presence of a thermodynamic 
singularity. In particular, optical absorption studies of the transition usually yield smoothly sigmoidal behavior [p), 
although this might be attributed to finite size effects or to the sequence heterogeneity of DNA. However, the PS 
models do not incorporate an additional important effect, namely elastic strain. As shown in Figure [j], because 
of the interwinding of the two strands of duplexed DNA a denaturation bubble induces extra elastic twist in the 
undenatured portions of the strand. For circular DNA, this strain energy cannot be relieved in the absence of special 

■ cutting enzymes, and progessive supercoiling can be expected to arrest the denaturing transition. The thermal 
denaturation of circular DNA is well-known to be sensitive to the degree of over- or underwinding of the loop before it 

■ is closed . In the biophysics literature, sigmoidal thermodynamic behavior is also connected with the development 
C5 of supercoiling ||] . 

It is the aim of this letter to demonstrate that when supercoling is included in the PS model, the thermodynamic 
singularity is, indeed, smeared out or seriously weakened. The starting point is the inclusion of distributed twist 
energy M in the PS Hamiltonian: 
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The first term in (Q) is the distributed twist elastic energy. The untwisting imposed by the unbinding of N2 base pairs 
is assumed to produce a twist proportional to N2 on the remaining N\ intact bases. Here, a quantifies the extent to 
which the circular chain is underwound or overwound. The limit a = 1 corresponds to complete unwinding of the 

O ■ DNA. When a < 0, the DNA duplex is overwound. The constant K is proportional to the torsional rigidity of double 
helical DNA (K ~ 10 -12 erg). The second and third terms describe, respectively, the base-pairing energy (eo ~ 2 

^ 1 kcal/mole) and the excluded volume interaction mentioned above. 
Is^j ' We will examine the phase behavior of this Generalized Poland Scheraga (GPS) model in the grand canonical 
ensemble, in which the partition function Z K '(z\, Z2) depends on the fugacities, z\ and Z2, of, respectively, intact and 

C$ ■ broken base pairs as: 
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Z K ( Zl ,z 2 )= J2 z i lz 2 2 Z^, N2 (2) 
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The partition function, Zj$ N ^ of a chain with prescribed numbers, N\ : 2 of respectively intact and broken base pairs 
is related to the corresponding partition function Z^ = ^ of the stress-free PS model by 
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Applying Cauchy's Theorem to Eq. (g), we obtain 

The functional form of Z K=0 (zi, z 2 ) near the thermodynamic singularity of the PS models adopts one of two different 
forms: 

Z K =°( Zl , z 2 ) oc ((z e (I) - z x ) + c(I)t + b(I)(z e (I) - z 2 ) - a(I)(z e (I) - z^y 1 (5a) 
Z K=0 ( Zl , z 2 ) oc ((z e {II) - z x ) + c(II)t + a(II){z e (II) - z 2 )P")- 1 (5b) 

where z e {I,II), c(I,II), a(I,II) and b(I,II) are positive numbers that depend on /3eo- The parameter t is the 
reduced temperature, which vanishes at the melting point. The exponents pi > 1 and < pu < 1 play a key role in 
the thermodynamic properties of the melting transition. 

Case I (Eq. (|5a|)) corresponds to PS models in which excluded volume interactions are included, both for the more 
flexible single strands and for double strands. These interactions lead to a first order melting transition. However, in 



case II (Eq. (5b)) excluded volume interactions are included only between single strands. In Case II, at the critical 
temperature t = a pole in the complex plane merges with the branch cut starting at z e {II) 1 corresponding to a 
continuous phase transition. The mean length of a denaturation bubble diverges at t = 0, while the correlation length, 
£(t), diverges as t 1 ^ 11 , and the specific heat exponent a = 2 — 1/pu. The power pu thus plays the role of a critical 
exponent, and hyperscaling is obeyed with the dimensionality d equal to one. In Case I, the mean size of a bubble 
remains finite at the temperature. The melting transition is first order. 
Performing the integration over z\ in Eq. (|J) we obtain 



Z« uN2 oc d> dzexp [-Nf UI (z, M)] (6) 

with M = N 2 /N the fraction of all pairs that are broken, and z = z 2 . The function fiji is the sum of entropic and 
enthalpic terms: 

f U i(z,M) =M\nz + (l-M)\n(zff I °(0,z)) + pK(M - a) 3 /(l- M) (7) 

In the thermodynamic limit N — > oo, but with M finite, the partition function is dominated by the minimum of fiji 
with respect to both z and M. Effectively fjjj(z, M) plays the role of a Landau variational free energy per base pair. 
This leads to the following coupled equations of state for M(t): 
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We first focus on the case a = 0. Using Eq. (5b) for Z A=0 (0, z), we obtain the results for M(t) shown in Fig. |[ For 
Case II there is no thermodynamic singularity at t = 0. The function M(t) has a smooth, sigmoidal behavior, even 
for very small values of K, and there is no specific heat anomaly. Physically, this is due to the continuous arrest of 
the denaturation by the build-up of strain energy as more and more bubbles open. The transition is never completed, 
and the "denatured" case consists of large bubbles connected by tightly twisted double-stranded sections. For Case 
I, the thermodynamic singularity is not completely removed, but the first order phase transition is transformed into 
a weak, third order transition || (see Figure ||) . 

We now turn to the the effects of non-zero values for the underwinding parameter a on the melting transition. In the 
case of a transition that is first order in the absence of strain-related effects (Case I), we find that, with one exception to 
be outlined below, underwinding or overwinding do not change the qualitative effect of strain energy on the transition. 
The "renormalization" of the order of the transition is the same, although the transition temperature will change 
with a. Figure ^ displays the dependence of the denatured ratio, N 2 /N as a function of reduced temperature for 
different values of a. Of note is the fact that underwinding can lower the transition temperature below its value for 
K = 0. Underwinding also enhances denaturation at low temperatures, even though the full transition is frustrated. 
This qualitative result is in general accord with recent experimental observations [ |10| ] . 

Figure |5| displays the transition temperature of under- and overwound DNA as a function of a. Note the minimum 
in the transition temperature at a — I. At the special value a = 1, which corresponds to a strand of duplexed 
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DNA that has been underwound to the point that the two individual strands are completely unwound, the transition 
temperature is lowered, and the first order transition to full denaturing is recovered. 

It thus seems clear that inclusion of denaturing-induced stresses in the PS model for circular DNA either removes 
the thermodynamic singularity (Case II) or weakens it to such an extent that one would encounter great difficulty in 
observing it experimentally (Case I). The one exception is the introduction of "complete" underwinding, corresponding 
to underwinding parameter a = 1. However, an important limitation on the present model is that we have assumed 
that the strains created as the result of the appearance of denaturation bubbles are taken up entirely in the form of 
twist. It is well-known that torsional stess on undenatured portions of DNA also produces writhing, e. g., formation 
of heterogeneous three-dimensional interwound plectonemic structures. It is likely that the effect of writhing will be 
to suppress the re-appearance of the first order transition for a = 1, because plectonemic structures are inherently 
inimical to strand separation. Such structures might account for the non-monotonicity of melting temperature as a 
function of a observed in 110(1. 
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FIG. 1. The creation of a denaturation bubble induces extra twist in the undenatured portions of DNA. 
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FIG. 2. Dependence of the fraction of broken pairs, M = N2/N on the reduced temperature, t, for Case II with pu = 0.5 
and cr = 0. The dashed curve shows M(t) for the K = case, i.e. with denaturation-induced twist energy fully relaxed. The 
solid curve shows M(t) for f3K = 1. In the second case there is no thermodynamic singularity. The melting transition in the 
strain-unaffected case is at t = 0. 
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FIG. 3. Dependence of the fraction of broken pairs, M = N2/N, on the reduced temperature, t, for Case I with pi = 1.5. 
The dashed curve shows M (t) for K — 0, and the solid curve shows M(t) for f3K = 1. In the second case, there is a third order 
thermodynamic instability. 
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FIG. 4. Dependence of the fraction of broken pairs, M = N2/N, on reduced temperature, t, for different values of the degree 
of over- and underwinding for Case I with pi — 1.5 and f3K — 1. For cr / 1 there is a third-order thermodynamic singularity. 
For a — 1 the first-order singularity of standard Case I PS models is recovered. The dashed curve is of the fraction of broken 
pairs when K = 0. 
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FIG. 5. A plot of the transition temperature, t(a), of under- and overwound DNA as a function of the degree, a, of over- or 
underwinding, for Case I wilth pi = 1.5 and f3K — 1. When a = 1, the interwinding of the two strands of duplexed DNA has 
been removed. 
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